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C — is an index indicating the latitude band 
O— (cap) indicates ocean 
s— indicates surface 
w— indicates water 
cloud fraction 

is the fraction of the latitude band that is land (L), ocean (O), land ice (LI) or 
ocean ice (OI) 

is the function which relates cloud optical thickness to cloud albedo A c 
arc through which the Earth rotates in time t 
the solar irradiance on day d 

the slope function relating surface albedo to top of the atmosphere 
the time 

time interval measured from noon to either sunrise or sunset 
is the earth’s rotation period (24 hours) with respect to the sun 
is weight function 
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Greek letters (all lower case) 


delta 5 is the solar declination 

theta 0 o is the solar zenith angle 

mu fi Q is the cosine of the solar zenith angle 

phi 0 is the latitude 



INTRODUCTION 


Satellite observations of the earth’s albedo are constrained due to the characteristics of the satellite 
orbits. Most of our presently existing albedo data have been collected from satellites situated in sun- 
synchronous, near polar orbits. Such orbits allow observations at a fixed time of day for each latitude 
observed. The satellite data are then used in forming monthly averaged values. The assumption usu- 
ally made in forming such averages is that the albedo observed must be modified if an average 
daily albedo is desired. 

It has been shown that the albedo of clear and cloudy atmospheres do depend on the solar zenith 
angle. In the present study, the angular models of Nack and Curran (1978) have been used. They 
parameterized the radiative transfer calculations of Dave and Braslau (1975). They also set up a rapid 
algorithm to transform surface albedo to surface-atmosphere albedo and irradiance and to obtain tem- 
poral averages. An updated version of their algorithm was used to compute corrected monthly aver- 
aged zonal and global albedos based on Nimbus ERB-7 Wide-Field-of-View albedo measurements. 

The corrected results are compared with both the original measurements and with results from an 
alternate correction scheme used by the ERB Nimbus experiment team and based on empirical direc- 
tional reflection data. 

In recent years analysts have commonly used empirical directional reflectance functions (Raschke 
et al., 1973, Taylor and Stowe, 1984, and Brooks and Minis, 1984) to transform the observed albedos 
into daily averaged albedos. The directional reflectance function for a given scene type relates the 
overhead sun albedo to the scene albedo for any other solar zenith angle. The averaging algorithm 
used by the Nimbus-7 ERB Experiment Team is discussed in Kyle and Vasanth (1986). The Nimbus-7 
ERB experiment and the data processing algorithms used are described in Jacobowitz et al. (1984). 
These empirical reflectance functions refer to the top of the atmosphere albedo and therefore bypass 
atmospheric transmission problems and calculations. At present however, there is some uncertainty 
both as to the accuracy of the various empirical reflectance functions and of how to best use them to 
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obtain an optimal daily averaged albedo from the measurements of sun-synchronous satellites. In this 
paper we will use the measured Nimbus-7 ERB Wide-Field-of-View orbital albedos to calculate 
monthly averaged, latitude band albedos by the updated method of Nack and Curran (1978) and 
Dhuria (1981). We will then compare our results to the monthly averaged albedos calculated by the 
ERB Experiment Team’s algorithm. 

COMPUTATION OF TOP OF THE ATMOSPHERE ALBEDOS 

A brief review of the Nack and Curran (1978) procedures and updated versions, Dhuria (1981), is 
presented here. Dave and Braslou (1974 and 1975) performed radiation transfer calculations for a 
realistic atmosphere with and without a stratus cloud layer. Nack and Curran (1978) parametrized 
their results and set up an algorithm to rapidly calculate the top of the atmosphere albedo for a cloud- 
free, a cloudy atmosphere, and a mixed partly cloudy /partly cloud-free atmosphere. The algorithm 
includes an isotropic surface reflectivity and can handle a cloud layer of arbitrary thickness. Temperal 
averaging was also included so that the average daytime albedo could be calculated. The algorithm 
does require as input the percentage of cloud cover. Curran et al. (1978) applied this algorithm to 
Nimbus-6 ERB latitude band monthly averaged observed albedos for the period July 1975 through 
April 1977. They first used the algorithm to estimate cloud fractions at each latitude from the 
observed albedos. Then using these cloud fractions the daily averaged albedos were calculated. 

The albedo at the top of the atmosphere can be expressed in terms of cloud albedo and cloud-free 
albedo as: 

A = A S (l-f c ) + A c f c (1) 

where 

A s = albedo at the top of cloud free atmosphere 

A c = albedo at the top of cloudy atmosphere 

f c = cloud fraction 
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Given measured values of A the cloud fraction can be obtained by first calculating A s and A c then 
solving for f c . To compute the daily averaged value of the albedo A d ^ at latitude 0 the time from 
sunrise until local noon (t d ^) is first computed from the equation of astronomy describing the apparent 
motion of the sun: 

0 o (t) = cos-Vo(t) 

ti 0 ( t) = sin6 d sin0 + cos6 d cos0 cos H(t) (2) 

H(t) = 27rt/T 

where 

0 Q (i) is solar zenith angle 

£ is an index indicating the latitude band 

0 is the latitude 

6 d is the solar declination on the dth day 

T is the earth’s rotation period (24 hours) with respect to the sun 

t is the time 

t d f time interval measured from noon to either sunrise or sunset for given S d and latitude 
band £. 

Averaged daily values of A s and A c at latitude £ are given by: 

\» - *»<*• <*> - /_'* A ‘ (a "'- ^ <3) 
= A c (a d „ t d , , r c ). ± 1'^ A t (a d „ 4 0 <t>, r c ) (4) 

where a df is the average surface albedo 
r c is cloud optical thickness 
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The surface albedo is assumed to be independent of the wave length and Lambertan reflection is 
adopted. The latitude band average, climatological surface albedo, a d ^, were calculated by using 
climatological snow and ice cover, and geographical ocean and land fractions and the relationship. 

aW>, month) = (f L1 + f OI )a r + (f L - f LI )a L (5) 

+ (1 _ fL _ foi) a w 

where fj is the fraction of the latitude band that is land (L), ocean (O), land ice (LI) or 
ocean ice (OI). 

a w is the albedo of water at local noon 

In these calculations only latitude band averages are treated. For this reason and for simplicity con- 
stant values of aj and a L were used: 
a, = 0.75 
a L = 0.18 

The fractions of land, ocean, land ice, and ocean ice were taken from tables 1, 2, and 3 of Curran 
et al (1978). Their table 8 also lists the ocean albedo, a w , as a function of latitude and time of year. 

At the equator a^^ = 0.06 but in the Arctic Ocean it has a value of 0.10 to 0.11 in the summer and 
of 0.23 to 0.25 in the early spring or late fall. Now jZ df is the average value of /* 0 ( t) at latitude (£) while 
(d) represents the day of the year: 

Nt = 57- t de 
2t d ‘ J-'dt 

Next to fit the results of Dave and Braslau (1974; 1975) linear functions of A g and A c are con- 
sidered as: 

A s (a, n 0 (t)) = m s (n Q ) a + A s (o, /i Q ) (6a) 

A c (a, n Q (t), t c ) = m c (ji Q ) g (r c ) a + A c (o, n 0 , t c ) (6b) 

m s and m c are slope functions and a is the surface albedo. 

Terms dependent on higher powers of a were neglected since the fit resulting from the linear function 
seemed satisfactory. 
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The slope factor for 6a depends on /i 0 ( t). The slope of 6b has been separated into a factor dependent 
on /i 0 ( t) and a factor g(r c ) dependent on the cloud optical thickness. 


In the present theory g(3.35) = 1. The factor g(r c ) adjusts the equation for the cloud optical 
thickness. It is defined by the relationship: 


g(r c ) = A c (0,0, 3.35)/ A c (0,0,r c ) 

A c (0,1, r c ) = 0.48 was chosen as the climatological tunning factor of global cloud albedo. To speed 
integration procedures, the slope and intercept functions of Eqn. 6 are expressed as functions of n Q in 
equations 7-10. The coefficiencies in the series expansion were obtained by a least square, minimum 
difference, fit to the calculations of Dave and Braslou. 


A s (0,^ o (t)) = 0.35057 - 1.0933jU Od£ + 1.6599ju^ 

- 1.1897 nl n + 0.32105 /£ n 
°d £ °d£ 

and A c (0, fi 0 (t), r c ) = 0.73 - 0.25^ 

m s d G* 0 (0) = 0.31876 + 1.2638 ^ - 1.368 ^ + 1.50833 ^ 
m c d (m 0 (0) = 0.16325 + 0.3633 ^ - 0.02501 ^ 


( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 


ALBEDO COMPUTATIONS 


Equations (3) and (4) now take the form: 

Asj/M’MoOdt)) = f dt t m s (m 0 + A s(0>Mo)] (M 0 (t)//i(t))dt (11a) 

Zl df J-t d( 


The monthly averaged latitude band albedos A. and A r are obtained from the daily albedos 

s ml ml 

weighted by the daily insolation or equivalently: 


£ 

ii 



where I d is the solar irradiance on day (d). Thus 



/ JL \ . 

JL \ 


A Sm , - ( £ A, *(,,) / | 

E w df 

\d=l 

| .02) 


/ n ^ 


A <W = ( E A c W d ,j / | 

E w df 

\d=l ) 

| (13) 


where n = the number of days in the month, and m indicates the month. 
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The average monthly latitudinal cloud fraction, required to complete the albedo A in equation (1) is 
obtained from: 


c m£ 


(A -S- A « ,/(A c* mS - A U' 


where A mf is the measured monthly averaged albedo. We used Nimbus-7 ERB values of A m ^ for the 


period January 1979-October 1983. The clear and cloudy sky albedos. A* , are the latitude band 
average values at the time of observation; this is close to local noon at most latitudes. 


ANALYSIS 


Three year-averaged measured orbital albedos, for 1979, 1980, 1981, are shown in Figure 1 in a 
latitude band averages versus time plot. Minimum albedos of about 22 percent occur in the tropics 
with steadily increasing albedos towards the poles. Maximum values of 70 to 80 percent occur in the 
polar regions in the spring and fall. The albedos shown during the polar night are found by interpola- 
tion and are in a sense fictitious, since the sun is below the horizon. The equivalent daily averaged 
albedos calculated by the ERB Experiment Team are shown in Figure 2 while our results are shown 
in Figure 3. The ratio of our averaged albedos to the orbital measured albedos is plotted in Figure 4, 
while Figure 5 shows the ratio of our results to those of the ERB Experiment Team. 

The daily averaged albedos, Figure 2 and Figure 3, show a pattern similar to Figure 1, however, the 
differences between the measured and daily averaged albedos are significant. Figure 6 is a graph of 
local time versus the latitude of the Nimbus-7 subsatellite point. Most of the Earth is observed near 
local noon (11:00 a. m. -1:00 p.m.). For a fixed scene, albedos are almost always a minimum at local 
noon and increase as the sun gets lower in the sky. Thus, daily averaged albedos should be larger 
than observed values at low and mid latitudes but smaller near the terminate which always lies in the 
polar regions for the Nimbus-7. 

This behavior is illustrated in Figure 4, where our computed daily averaged albedo is expressed as 
percent of the measured albedo. We calculate the daily average albedo to be between 110 percent and 
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114 percent of the measured, noon time, albedo in the tropics but it drops to between 94 percent and 
98 percent of the polar sunrise or sunset albedos. Taylor and Stowe’s (1984) data indicates that the 
noon time-terminate albedo range in the tropics should be a factor of 2 or larger. However, the aver- 
aged albedos are insolated weighted and the most intense irradiance occurs near local noon. In the 
polar regions the sun is always low in the sky. Thus, the measured albedo at 6 a.m. or 6 p.m. will 
not vary too much from that measured at local noon. In the tropics however, the solar zenith angle 
will vary from between (0-30°) at local noon to 90° at sunrise or sunset. Correspondingly for a fixed 
scene the measured albedo will increase from local noon to sunset by between 50 percent and 200 
percent depending on the scene. 

There is no established optimal procedure for estimating the daily average albedo from a once a day 
measurement. The daily averaging algorithm used by the ERB Nimbus Experiment Team (Kyle and 
Vasanth, 1986), is known to be simplistic but tests indicate that it is reasonably accurate. Its two chief 
defects are that it uses a global, scene independent algorithm to describe the variation of albedo with 
solar zenith angle and that it is based on Nimbus-2 and 3 MRIR albedo measurements which are not 
as accurate as more recent satellite albedo measurement. 

Our computed daily averages are compared to the Experiment Team’s daily average estimates in 
Figure 5. Our values are expressed as percentages of the Team’s values. The major differences 
between the two daily averaging algorithms should be kept in mind: 

• The ERB Team algorithm is scene independent and varies only with solar declination and time of 
observation. Our algorithm, however, is dependent on the zonally averaged scene (percent of land, 
ocean, ice and cloud in the latitude band). 

• The ERB Team algorithm starts with the measured albedo and simply applies the daily averaging 
factor. Our algorithm uses climatological surface and cloud albedos combined with radiation 
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transfer theory to obtain the top of the atmosphere albedo. The measured albedos are used to 
estimate the amount of cloud present. 

In the tropics and mid latitudes the agreement is generally good (+2 percent). However, our values 
drop 6 percent below the Team values in a number of places. These regions are associated with 
minimum cloudiness as indicated in Figure 7. This indicates a definite defect in the Team algorithm 
which does not adjust for variations in the scene. 

In the polar regions the difference between the two algorithms range up to 10 percent but this should 
be expected since neither algorithm treats these regions properly. The directional reflectance for snow 
and ice is considerably different from the scene independent model used by the ERB Team. On the 
other hand our algorithm does a poor job of estimating cloud amounts in the polar region. Its zonally 
averaged correction factors are therefore less reliable in these regions. 

The climatological values used in our model assumes mean global clear sky land, ocean and snow 
albedos and fixed overcast albedos. Although our daily averaged albedo values appear reasonable in 
all zones some values of the derived cloud cover in the polar regions appear too low. For instance, a 
value of 14 percent cloud cover at 60° N latitude in April is much too low. The problem seems to 
arise from the use of the mean global scene albedos. Use of seasonal and zonal scene albedos would 
markedly improve our cloud values in regions like this. The low mean sun angles in the polar regions 
keeps our daily averaged albedo from straying too far from the true value. 

CONCLUSIONS 

The use of a global, scene independent directional reflectance function to calculate daily averaged 
albedos does not yield accurate seasonal and zonal variations. The problem is particularly bothersome 
in the tropics and mid latitudes where cloudy and clear zonal bands occur. The positions and inten- 
sities of these bands also vary seasonally. 
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A heuristic algorithm based on radiative transfer theory has been used to study the problem. The daily 
averaging functioning accounts for both zonal geographical differences and for varying cloud cover. It 
has been applied to three years of Nimbus-7 ERB Wide-Field-of-View measured albedo data. 

When compared to the global, scene independent algorithm used by the Nimbus ERB Experiment 
Team it yields more realistic results in the tropics and mid latitudes. In the polar regions both 
algorithms have noticeable faults. 

Based on our results we recommend that an improved daily averaging algorithm be set up to calculate 
daily averaged albedos and Net radiation from the Nimbus-7 ERB MATRIX tapes. These tapes con- 
tain both the observed and daily averaged albedos. The algorithm described in this paper would be 
useful for this purpose. However, first it should be modified to use seasonal and zonal climatological 
surface and cloud albedos. 
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Figure 1. Nimbus-7 ERB Wide-Field-of-View measured climatological latitude band average albedos. 
Data covering three years, January 1979-December 1981, have been combined to form an average 
climatological year. Latitude band averages are plotted against month of the year. The albedos are 

shown in percent. 


11 



* 

* 

# 

L 

A 

T 

T 

I 

T 

U 

D 

E 

x 

x 

X 


90.00 r 

85.50- 

81.00- 

76.50- 

72.00 - 

67.50 - 

63.00 6 

58.50 = 

54.00- 

49.50 = 

45.00 - 

40.50 = 

36.00 - 

31.50 = 

27.00 = 

22.50 - 

10.00 - 

13.50 - 

9.00- 

4.50- 

0.00 - 

-4.50 - 

-9.00 - 
-13.50 - 
-18.00- 
-22.50- 
-27.00- 
-31.50- 
-36.00 = 
-40.50 - 
-45.00 - 
-49.50 - 
-54.00 z 
-58.50- 
-63.00 f 
-67.50 - 
-72.00 S 
-76.50 =- 
-81.00- 
-85.50 - 
-90.00* — 

JA^J 


Figure 2. C 
Experiment ' 



estimate the 




* 

* 

* 

L 

A 

T 

T 

I 

T 

U 

D 

E 

* 

* 

* 


90. DO 

85.50 

81.00 

76.50 

72.00 

67.50 


63.00 

58.50 

54.00 

49.50 

45.00 

40.50 

36.00 

31.50 

27.00 

22.50 

18.00 

13.50 

9.00 

4.50 

0.00 

-4.50 

-9.00 

-13.50 

-18.00 

-22.50 

-27.00 

-31.50 

-36.00 

-40.50 

-45.00 

-49.50 

-54.00 

-58.50 

-63.00 

-67.50 

-72.00 

-76.50 

-81.00 

-85.50 

-90.00 



jaKi FE?B MA& aA JlA juL Alft s$> o£9 n<W 


Figure 4. Ratio in percent of the calculated monthly averaged albedos (Fig. 3) and the measured 
albedos (Fig. 1) given in percent. Latitude band ratios are plotted against month of the year. 
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Figure 5. Ratio in percent of the calculated monthly averaged albedos (Fig. 3) to the monthly aver- 
aged albedos calculated by the ERB Experiment Team (Fig. 2). Latitude band ratios are plotted 

against month of the year. 
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Figure 6. Nimbus-7 local sampling times for nadir point as a function of latitude. 
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